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In this paper we present a comprehensive study of inclusive hadron production in DIS at 
low X. Properties of the hadron spectrum are different in different kinematic regions formed 
by three relevant momentum scales: photon virtuality Q^, hadron transverse momentum kx 



p^ and the saturation momentum Qa{x). We investigate each kinematic region and derive the 

H corresponding asymptotic formulas for the cross section at the leading logarithmic order. 

^I We also analyze the next-leading-order (NLO) corrections to the BFKL kernel that are 

T-H responsible for the momentum conservation. In particular, we establish the asymptotic 

^ behavior of the forward elastic dipole-nucleus scattering amplitude at high energies deeply 

f-^ in the saturation regime and a modification of the pomeron intercept. We study the nuclear 

fS , effect on the inclusive cross section using the nuclear modification factor and its logarithmic 

(-H derivative. We argue that the later is proportional to the difference between the anomalous 

dimension of the gluon distribution in nucleus and in proton and thus is a direct measure 

►^ of the coherence effects. To augment our arguments and present quantitative results we 

[£j performed numerical calculations in the kinematic region that may be accessible by the 

r future DIS experiments. 
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I. INTRODUCTION 

In the last decade we have learned a great deal about gluon saturation/color glass condensate 
[TVI27J thanks to the relativistic dAu and AuAu program at RHIC. The future DIS programs at EIC 
and LHeC promise to provide even more detailed information about structure of the nuclear matter 
at low X. How successful that program will be depends a lot on our ability to pinpoint the processes 
that are most sensitive to the low-x regime. In this paper we study one such process - inclusive 
hadron production in eA scattering. It has been a subject of intense theoretical investigation 
over the past decade [28H37J and has proved to be a powerful tool in dA collisions at RHIC. 
On the one hand, we expect that p{d)A and eA processes have very much in common due to 
the Pomerantchuk theorem, that states that all high energy scattering processes are mediated by 
exchange of a collective gluon state - known as pomeron - that has vacuum quantum numbers. On 
the other hand, proton wave function is characterized by a soft, non-perturbative scale, whereas the 



virtual photon wave function can be calculated using the perturbation theory and is characterized 
by virtuality Q^ . A possibility to dial Q^ is a great advantage of DIS. Our main goal in this paper 
is to provide a thorough analysis of the inclusive hadron production in various kinematic regions 
characterized by three dimensional scales: photon virtuality Q^, hadron momentum k^ and the 
saturation momentum Qs and to produce numerical predictions for both novel and well-known 
quantities that can be tested at EIC and/or LHeC. 

Our paper is organized as follows. In Sec. |ll] we use the dipole model [38] to relate the DIS 
^* A cross section to that of the color dipole qq + A. The ^* A differential cross section can be 
expressed in a factorized form as a product of the light-cone wave function of the virtual photon 



7* and qq^ A differential cross section. In Sec. Ill we review the properties of the BFKL pomeron 



Sol HO] and the unintegrated gluon distribution function at LO, particularly we emphasize the 



leading logarithmic asymptotics. These are used in Sec. |IV] to derive the asymptotic properties of 
gluon production in dipole-nucleus scattering in various kinematic regions. In Sec. |V] the result is 
further generalized to the case of LO gluon production in DIS. 

The NLO corrections to the inclusive hadron production are rather complex. These include NLO 
correction to the BFKL kernel [4H H2] , |43fl52j , running coupling corrections |53H60j and energy 



conservation |61H63j corrections to BK |21H24j . It has been argued in [63j that energy conservation 



is the most important phenomenological effect beyond the LO. Therefore, in Sec. VI we investigate 
the role of this effect on inclusive hadron production. In our calculations we rely on a phenomeno- 
logical approach suggested in |641 165] where a modified BK (mBK) equation that satisfies energy 
conservation was derived. It was utilized in |61t [66] to calculate the NLO corrections to the total 
DIS cross section. niBK equation serves as the basis for our NLO calculations. First, we derive 
the dipole scattering amplitude in dilute and saturation regimes; the corresponding expressions are 



given by (70) and (83) respectively. We argue that the energy conservation effects decrease the 
energy dependence of the saturation momentum. These results are used for computation of dipole 
density in various asymptotic regimes. Similarly to our analysis of LO case, we explore the NLO 
gluon production first for dipole — nucleus process and then for DIS scattering. 

It is very instructive to know how the DIS on a heavy nucleus is different from DIS on a proton 
at low X. Had the coherence length been short, of the order of the proton radius, the hadron 
production in ^* A would have been equal the incoherent sum oi A 'y*N processes. However, since 
the coherence length is larger than the nuclear radius, the entire process is coherent. Because it is 
interesting to compare the coherent and incoherent regimes, one introduces the nuclear modification 
factor (NMF) R that calibrates the cross section in 'y*A with that of 7*A^ rescaled by atomic weight 



A. Sec. |VII| is devoted to the study of the properties of this quantity as a function of the hadron 
transverse momentum, photon virtuality and atomic weight. 

We expect that at EIC/LHeC kinematic region the low- a; evolution effects start to play an 
important role rendering the anomalous dimensions dependent on atomic weight. This manifests 
itself in inclusive hadron production in dA collisions at RHIC as the transition from the Cronin 
enhancement at mid-rapidity to suppression of the NMF at forward rapidities even at kx > Qs- In 
order to evaluate how steep is the dependence of the NMF on rapidity, we introduce a new observ- 



able J, defined as the logarithmic derivative of R, viz. dlnR/dy. We demonstrate in Sec. VII that 
at kx ^ Qs, J is proportional to the difference of the anomalous dimensions of the gluon distribu- 
tion in nucleus and in proton. Without the low-x evolution one expect J to vanish. However, due 
to the low-x evolution J acquires a finite negative value. Therefore, J can serve as a direct probe 
of the effect of the slow- a; evolution on the nuclear gluon distribution function. 



The numerical computations are presented in Sec. VIII We use the bCGC model [67] for the 
dipole-nucleus forward scattering amplitude, albeit with the simplified 6-dependence. In Fig. [5] 
we plot d'^ F2/d In k'^dy as a function of photon virtuality Q^ and hadron transverse momentum 
kx and rapidity y = ln(l/xp).* In order to emphasize the role played by the NLO effects we 
exhibit both LO and NLO results in each plot for the structure function. In Fig. [5] we see that 
the NLO calculation yields much smaller cross section for inclusive hadron production than the 
LO one. Additionally, its functional dependence on kx, Q^ and y is substantially weaker in NLO 
than in LO. This is in accordance with our observation in Sec. IVII that NLO correction reduces 
the anomalous dimension of the gluon distribution. Interestingly, most of the NLO effect cancels 
in the NMF which appears to be a robust quantity in this respect. This indicates that the energy 
conservation effect factors out to a large extent from the inclusive cross section. 

The NMF shown in Fig. [7] displays a number of interesting features. First, the NMF is strongly 
suppressed at small /ct's but exhibits an enhancement toward higher fc^'s where the Cronin effect 
(R > 1) is observed. This seems to be in contrast with pA collisions [34J where the Cronin effect 
gives way to suppression of NMF at all kx's as the hadron rapidity increases. This is the result 
of the linear evolution in the rapidity interval between the virtual photon and the hadron. This 
evolution produces dipoles of different sizes that scatter in the nucleus with different amplitudes. 
At small kx large dipoles, on which the gluon saturation effects are stronger, dominate the cross 



* We use the a;ip notation borrowed from the diffractive DIS where it denotes the momentum fraction carried by the 
pomeron. It does not have this simple interpretation in our case because the interaction is inelastic. 



section, whereas at higher kT smaller dipoles contribute to the NMF enhancement. Second, we 
observe a relatively weak A-dependence. This is also a result of the averaging over different dipoles. 



Third, we note a peculiar Q dependence that is explained in Sec. VIII 

To investigate the rapidity dependence in more detail we plot the logarithmic slope of the nuclear 
modification factor J on Fig. M (for dipole-nucleus scattering). We see that it is negative for the 
entire kinematic region indicating the graduate suppression of the NMF towards large rapidities. 



This is in agreement with our arguments in Sec. VII We argue that J is directly proportional to 



the difference between the anomalous dimensions of the gluon distribution function in the nucleus 
and in proton. Hence we believe that measuring J is a great tool for exploring the low-x regime of 
QCD. 

We summarize our results in Sec. IIXI 

II. FROM -f*A TO qq + A SCATTERING 

The dominant contribution to the inclusive hadron production in DIS at low-x, at rapidities 
away from the virtual photon and nucleus fragmentation regions, comes from the fragmentation 
of fast s-channel gluons [2j. The cross section for inclusive production of a gluon of transverse 
momentum k at rapidity y in deep inelastic scattering can be represented as an integral in the 
configuration space [68]:^ 



(Pkdy 27r2 J Jo ' ' d?kdy 

where the virtual photon wave function $ describes splitting of a photon of virtuality Q^ into qq 
color dipole. It is given by 

$(r,z)(Q,r,z) = |*^*(Q,r,z)|2 + |xI/f(Q,r,z)|2, (2a) 

/ 
$T(r,z)(Q,r,z) = 2N,y^^^{a'Kl{Ta)[z'' + (1 - zf\ + m)Kl{Ta)} , (2b) 

/ 

$L(r, z){Q, r, z) = 2N, V "^AQ^z'{l - zfKl{ra) . (2c) 

^ TT 

Here a? = Q'^z{l — z) + ml, a'^^ = e^2;?/(47r), with Zf being electric charge of quark / in the units 
of electron charge e. The cross section for inclusive gluon production in dipole-nucleus scattering 



^ We use the notation k^ — k'^ = k^, where k is a vector transverse to the coUision axis. 



reads [29] 



j dH I dV'e-^-E' [VI,Ng{i', b', y)] [V;Mh I, Y - y)] , (3) 



(Pkdy vr^ /c^ 

Here the dipole density n{r, r' ,Y — y)d'^r' is the number of daughter dipoles of size r' in the interval 
d'^r' produced by a parent dipole of size r at the relative rapidity Y — y |4h6J. It satisfies the BFKL 
equation [39^ HQ] with the initial condition 

n(r,r',0) =(5(r-r'). (4) 

At the leading logarithmic order, the corresponding solution is |39| HO] 



oo 



with the eigevalue function x given by 

X(z^) = ^(1) - ^^ Q - ii^) - ^V Q + ^^^) , (6) 

where ip{iy) = r'{u)/r{u). 

Let /(r,r',y) be the particular solution of the two-dimensional Poisson equation 

V2,/(r,r',y)=n(r,r',y). (7) 

Employing ([5]) we derive the Melin representation of / 

It is convenient to write Q as a convolution in the momentum space. To this end we introduce 
the Fourier-image of / with respect to the second argument: 

/(£,q,.)= /dV'e-'9-f /(t,t'.„) ^ ^ r d.."'-v<-)» (r^)"- '; ^\: '■'',, (9) 

and the unintegrated gluon distribution function of the nucleus [21 [29] 

V'Ads, y) = ^-^ j dH j d\ e-'^-'- V^Ncir, b, y) . (10) 

Ng{l, b, y) is the forward scattering amplitude of a color gluon (or adjoint) dipole r on the nucleus 
at impact parameter b at the relative rapidity y. It obeys the BK equation [211 [23] and its properties 
are discussed in the next section. Using ^ and ( 10 ) in ^ we get 



dai'i+^{k,y;r) _ Aa^ 
d^kdy irk^ 



ld^p^A{p,y)f{r,k-p,Y-y). (11) 



III. LOGARITHMIC APPROXIMATIONS 



A. Asymptotic expressions for / 



It is worthwhile to hst here the asymptotic formulas for / in various kinematic regions (we 
follow notations of |69H71| were more details can be found). 



1. asH ^ In ^. In this case the eigenfunction mh can be expanded near its minimum x 
2 In 2 — 7(^(3)i^^. Expression under the z^-integral in (jsl) has a saddle point at 



= H2M . (12) 

''' 14C(3)Q.y ^ ' 



In this approximation integration over z^ in ([8| produces 



in^ra 



/(r,q,y) = - , ^ \^^ e^-'^' -^)y e' ^m^v , (13) 

with aP - 1 = 4asln2. 

2. rq < 2 and In — ^ a^y. In this region, the leading contribution to the z^-integral stems from 
the pole at iu = 1/2. Approximating the eigenfunction as x ^ 1/(1 ~ 2zi/) and employing 
the saddle point method in ^ again yields 

The saddle point is 



2-.P = 1 - Jff ■ (15) 

3. rq > 2 and In ^ ^ asH- Now, another pole in x dominates, x ~ 1/(1 + 2iz^) with the result 
for / 

/fcq,y) = ^- L^e^V^^I^l^ (16) 

and for the saddle point 



2^z^sp = -l + J^. (17) 



B. Properties of ipA 



Unintegrated gluon distribution cpA is defined by (10). NG{l,y,h) stands for tlie forward elastic 
gluon dipole scattering amplitude. At large Nc, the gluon dipole is equivalent to two qq dipoles 
each of which scatters with amplitude A^(r, y, b). Therefore, 

Neil, b, y) = 2N{v, b, y) - N\i, b, y) (18) 

The qq scattering amplitude satisfies the BK equation |2H EH] and its properties are well-known. 
Initial condition for the BK equation is the Glauber-Mueller formula |38j for the forward scattering 
amplitude N oi a qq color dipole on the nucleus: 

iV(r,b,0) = l-e-5E''3?o. (19) 

The gluon saturation momentum [1] at initial rapidity y = 0, which corresponds to the Bjorken 
variable xq such that y = In — , is related to gluon distribution function xG at x = xq as 

Qlo = ^^?^Pn\^) ^oG(xo, l/r^) , (20) 

where p is the nuclear density, T(b) is the nuclear thickness function as a function of the impact 
parameter b. The gluon distribution function at the leading order in Os, i.e. in the two-gluon 
exchange approximation, reads 

xG(x,l/r2) = ^^ln^, (21) 

with A being some non-perturbative momentum scale characterizing the nucleon's wave function. 



Using (19) in (18) we derive the initial condition for the gluon dipole scattering amplitude 

A^G(r,b,0) = l-e-3E''3?o_ (22) 

Let us now list some properties of the amplitude Ng, see [Ml EH] for details. 

1. At r ^ 1/QsO the BK equation reduces to the BFKL equation, which must be solved with 
the initial condition A^(r, b, 0) ss t^QIq/A. Small dipoles scatter independently, perforce 
Nq ~ 2A^. Thus, in this region 

iVG(r,b,y) = r dz^e^-xH^ {rQ,,f^^^'^ 1 1 + (1 - 2i^) In gf _ 
J-oo 8^ (1 - 2^l^)^ 

2. In particular, if r ^ l/<5so and In -4— » Usy the solution is 



rWsO 



3. For r ^ 1/QsO and asV ^ In -7^— we have 



NG{l,h,y) 



rQsO In 



WsO 

A 



4 v'l4C(3)7ra,y 






(25) 



4. The saturation region is characterized by the saturation momentum Qs{y)- With the double 
logarithmic accuracy it reads |72H74] 



Qsiy) = Qsoe'^^^y 

In the saturation region r > 1/Qs, solution to the BK equation is [72 

Af(r,b,y) = l-5oe-ii°'(^'^'), 



(26) 



(27) 



where Sq is a constant that can be determined by matching N from (27) with that of (23) 
at r = 2/Qsiy). Consequently, 



iVG(r,b,y) = l-SoVi'^'(^«=), 



(28) 



where we utilized (18). 



Eqs. (23)-(28) are derived with the logarithmic accuracy. We can calculate 99^ given by (10) in 



the same approximation as 



^A(k,y) 



Cf 

as(2vr)2 



i/fc d f d 



d'b I dr— (r — iVG(r,b,y) 



Cf f,2^ d 



as{2'iTyk J dr 



d'b—NG{r/k,h,y). 



(29) 



We stress that this formula holds only in the asymptotic regions specified in 1-4 above; still this 
is a very useful approximation as it captures the most essential features of the unintegrated gluon 
distribution. 



It is evident from (29), that in place of function A'^G'(r, b,y) it is convenient to use 
function NG{k,b,y) = NG{i^/k,h,y), where r = r/r. In particular, dNG{i'/k,b,y)/dr = 



-k'^dNG{k,h,y)/dk} Plugging (29) into (11) we obtain 



da55+^(k,y;r) _ asCp /" ,2, /" ^2 dNGip,b,y) ~ 



d?kdy TT^/c^ 



d'b / d'p 



d\n.{l/p) 



/(r, p - k, y - y) . 



(30) 



* We assumed in ( 20 1 that the fe-dependence factors out in the initial condition; perforce it factors out in the solution 



for heavy nuclei. Therefore, scattering amplitudes depend only on the absolute value of vector b. 



IV. PROPERTIES OF THE DIPOLE NUCLEUS CROSS SECTION 



To calculate the cross section for gluon production in dipole-nucleus scattering we need to 



evaluate the integral over the transverse momentum p in the right-hand-side of ( 30 ) . It convenient 
to consider the inclusive cross section at a fixed impact parameter b: 



5(k,y,b;r) 



da95+^(k,y;r) f UsCp 



d p an.,/-, /„„N fil,P-l,Y -y)- (31) 



d?kdyd?h \Tr'^k^ J J dln{l/p) 

When taking the p-integral with the logarithmic accuracy in various kinematic regions it is useful 



to keep in mind that (l28|),(l23|) imply that dNc/dln^l/p) ~ ln{Qs/p)exp{-ln^{Qs/p)} ifp<^Qs 
and dNG/dln{l/p) ~ Q'^Jp'^ if p > Qs, while (|l4|,(|l6| indicate that / ~ l/Zc^ if A; > 1/r and 
/~r2, if A;< 1/r. 



1. k ^ Qs ^ 2/r. Due to the strong ordering of the relevant scales we have 



(32) 



Using ( 23 ) we derive 



q/^^ ain(i/p) -' y_^"^ V ^ 

(ln#:)^/^ 



l+2iv 



Qso 



QsoY^"'" 1 1 + (1 - 2iz^) In ^ 1 + 2it/ 
8^ (l-2ii/)2 l-2w 

(33) 



^Tifj^*(-\/i^'"^)'^^''^ 



Thus, it follows upon substitution of (16) and ( |33[ ) into (32) and then into (30) that 



dai^+^{k,y;T) asCp 



cPkdy 



SvrSfc^ 



d^bQ 



sO 



(2a,y)5/4(2a,(y-y)lnf)i/4 



1 + W^ln%^ 



X e 



2^2MY-y)\n!f 2^2a,yln^ 



(34) 



2. /c » 2/r » Qs- Repeating the by now familiar procedure yields 

^ , /■' dNG(p,b,y) - 



We observe that the cross section in this case is exactly the same as (34). 



(35) 



3. Q, > A; > 2/r: 



--//*^W^/'''^'''-' 



With the help of (28) and (16) we get 



g = 2rS\ 



dv- 



1 



l + 2iu 



2as(Y-y) 
g l + 2iv 



(36) 
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Now, using r = In ^ in place of p 
Jk P 



HQs/k) 



drr e 



-T^-T{l + 2ilJ,) 



Q 



2ifj.+l 







) 1 



(38) 



Putting everything together yields 



d'^kdy 
4. Qs > 2/r > k: 



7r5/2A;2 



2j2a4r-y)ln=i 



(hi^)i/4(2a,(y-y))i/4 



(39) 



_ . , ' dNG{p,b,y) ~ 
g^lTT dpp f{T, p,Y -y) 

1 2/r dln{l/p) 



(40) 



This case is similar to the previous one except the the lower limit of the integral in (37), k, 



is now replaced by 1/r. However, for very large Qs, the integral over p is independent of the 



lower limit of integration as is clear from ( 38 ) . We conclude thereby that the cross section 



in this case coincides with ( 39 ) 



5. 2/r > /c > Qs-. 

^„ /■' dNG{p,b,y) ;, , ^ . ^ /•2/'- dNcipAy) ;, ^ , ,,,, 
g^lTTJ dpp ^,_^,., ,^_, f{i,k,Y -y)+27r dpp f{T,p,Y -y) (41) 



Qs dln{l/p) 



dln{l/p) 



The first of these integrals reads using (33) and (14) 



-''"-Hi^/'''^'^--) 



1 I / 2asy 1 Qso 

"32 (2a,y)5/4 / _ ,^^ ,, ^ ^^^ 



2 2 2^2a., In ^ 2^2as(Y-y) In ^ 



1/4 QsO'' ^ 



(42) 



The second one is done by substituting (24) and the integral form ([9]) (it is useful to note 
that dNG/dh\[l/p) ~ ^Nq) and then integrating over p in the leading log approximation 
(i.e. treating logp as a constant) followed by the saddle point integral over v. We have 






(i-^)^/^an^)^/Mi + 



=Q'o^ 



2as3/ i„ QsO 



In- 



In- 



A 



2(2a,y)3/4(2a,(y-y))3/4 



3V2^=^^°Qt<^eV2"-(^-^)i°^ (43) 
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Substitution of ( 42 ) and ( 43 ) into ( 30 ) gives for the cross section 



da--+^(k,,;r)_a.C^ / d^feQ^o-^ 



(Pkdy 



TT^fc^ 



(ln^)V4(in|,)i/4(i + ./,2a^ln2- 



QaO 



X e 



2y2^^2^1^^p2y2M>^^yy^ 



2(2a,y)3/4(2a,(y-y))3/4 

(ln^)V^(2a,(y-j,))V^ 
(2a.y)i/2(ln|:ln-^)V4 



1 + 



fcr rQsO - 



(44) 



6. 2/r>g, >A;: 






Repeating the steps leading to (43) and noting (26) we finally get 



da?«"+^(k,y;r) a.Cp 



(Pkdy 



TT^k'^ 



d%Qiy 






(45) 



25/2(2a,y)3/4(2a,(y-y))3/4 



,4V2a.«V2a, (>--?/) In Q^ 



(46) 



Eqs. (34)-(46) represent the dipole-nucleus inclusive cross section in all kinematic regions. 



V. GLUON PRODUCTION AT THE LEADING ORDER IN ASYMPTOTIC REGIONS 



The DIS inclusive cross section is obtained from the dipole-nucleus one using M. Integration 
over the dipole size r and momentum fraction z can be carried out for Q ^ A,m. In this case the 
largest contribution stems from the transversely polarized virtual photon. Setting nif = in ([2]) 
we write (II| as 

' ' f 



d^kdy 



^T^^Jd^rf dz Q'z{l - z) Kl (vQ^/^l 



['^-•-^fl^^^^H^ 



(47) 



At large Q the dominant contribution to the z-integral arises from z — )■ 0, 1. This corresponds 
to either quark or antiquark carrying most of the photon's energy. These limits are symmetric, 
therefore we can calculate the z-integral for z — )■ and multiply the result by 2. Thus, 



da^*^(k,y;Q) iY^Q^ 2a 



'OO 

em / , 2 



d'^kdy 



8Nr, 2a,ru f^ dr 2 dcT««+^ (k,y;r) 



vr^ 



3vr2Q2 3 7^/Q, ^4 



d^kdy 



(48) 
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where we took into account only three light quarks. To set the low limit of integration in (48) we 



noted that integrand in (47) peaks at rQ ~ 1/ ^Jz{l — z) >2. Upon substitution of (30) into (48) 
we get 



d'^kdy 97r5 Q2^2 J ^ " J^^q2 r^ J ^ '' ain(l/p) 

To determine the cross section for gluon production in DIS it is convenient to do integral over 
r before we integrate over u in /. We thus define an auxiliary function 



d{Q,p,y) 



4 r 



-fii,p,y) 



(50) 



Employing ^ in (50) we obtain the Mellin representation of d 



d{Q,p,y) = 77^ / due 



(51) 



,Qj {l-iiy)T{l + iu){2iu + iy 
Inasmuch as we are interested only in asymptotic behavior of d, which we will derive using the 
saddle-point approximation, we can write in view of ^ 



diQ,p,y) = ^f{2/Q,p,y)^ 



2 ^^sp 



(52) 



where fgp is a saddle point given by one of the formulas (12), (15), (17). In particular, using (13), 



(14) and (16) in (52) yields 



d{Q,p,y) 



d{Q,p,y) 



d{Q,q,y) 



Q 1 

4py/l47rC(3)a,y 






2P_ 

Q 



U/2 



A^/TT 



(^^qIo''^'' _^2^2asylnS. 



(2a,y)3/4(ln: 



1/4- 



Q' 



'^-P' ^2a.yln^'''' 



^2^2asyln^ 



Q'>P 
Q<.P 



(53) 

(54) 

(55) 



Inspecting (49), (50), (52), (30) and (31) we get 



da^'^{k,y-Q) ^ 4N,a,^Kdai^+'^{k,y;2/Q) 
d^kdy Dvr^ d'^kdy 

where we denoted by k the logarithmic (or constant) factor (1/2 — ivsp)~^. Explicitly, 

^^ max{fc,Q} \ 1/2 

K = 2\ ^_7y'^\ ,ifA;,Q»Q,; k = I , ii k,Q <^ Qs 

2as{Y-y) 



(56) 



(57) 



Eq. ( |56[ ) together with the expressions of the inclusive dipole-nucleus cross section derived in 
Sec. |V] provide the cross section for the inclusive gluon production in DIS at the leading logarithmic 
approximation. 



13 



VI. NLO BFKL EFFECTS: ENERGY CONSERVATION 
A. Dipole scattering amplitude 

As explained in the Introduction, one of the most important NLO effects is the momentum 
conservation. BK equation modified to account for the energy conservation reads [641 I65j 

'-^^^ = l(^ - 1) / ^'^'?^(^ ^^^^' ^' y^ + ^(^ - ^' ^' ^) + ^(^' ^' y^ 

-iV(r',b,y)iV(r-r',b,2/)}. (58) 
In this section we discuss solution to this equation in dilute and saturation regimes. 

1. Dilute regime 

Consider first the dilute regime. It is advantageous to represent N as the double Mellin trans- 
form 

iV(r, b, y) = / f^ / ^AA(7, b, c) \_ . , (59) 

where we introduced a new dimensionless variable ^ = ln(l/r^(5gg). The anomalous dimension 7 
is related to the Mellin variable v that we have used so far as 7 = 1/2 — iv, so that the LO BFKL 
eigenvalue function is x{^) = xiKl ~ l/^))) see Q. xiil^^) denotes the NLO BFKL eigenvalue 



function. In the dilute regime the A^^ term in the r.h.s. of (58) can be neglected. Substituting (59) 



into (58) one arrives at the following relation between the Mellin variables 

L^ = 2a,xi(7,w) = 2as(l-u;)x(i(7-l/2)) , (60) 

with the explicit solution for cj 

u = -^ . (61) 

This solution is plotted in Fig. [!} w diverges aX v = v* satisfying 2asx('^*) = ~1- As a^ — ^ 0, cj 
approaches the LO expression while v* — >■ ±00. At 7 — >■ 0, i.e. iu — )■ 1/2, x ~ 1/(1 — 2ii/) = I/27 



and (61) yields 



7(0;) = a, { ^ - 1 I . (62) 
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(a) 



^ 


a,=0.2 


""T 





(fc) 



FIG. 1: w(i/) for (a) a^ = 0.3 and (b) ois = 0.2. LO and NLO are represented by dashed (red) and solid 
(blue) lines respectively. Notice the different v ranges of the two plots. 



This can be used as a model of anomalous dimension that takes into account the energy conservation 
as suggested in [75l [76] . ^ 



Integrating ( 59 ) over uj we obtain 

iV(r,b,y) 



dv C^ e'^('^)^+^«-« 



(65) 



with w(z/) given by (61). Remembering that in the dilute regime (and Nc ^ 1) Nq = 2N, see (18), 



and using the same initial condition as in ( |23| ) we get 



(66) 



,l + 2a,xMJ "■''" Svr {1 - lii^y 

This integral can be taken in the double- logarithmic approximation (DLA), which corresponds to 
keeping only one of the poles of X) namely x(^) = 1/(1 ~ 2iz^). Denote 



Then, in the DLA 



<PiC,y) 



a.o 



7 + as 



HC, y) = 7—773 — ^ - 1/2 + iv)i ■ 
1 + 2asX(i^) 



(67) 



+ 7^ - e = 2v/^;^ - e(i + as) + - (7 - 70) 



2^3/2 
(^^^ 



(68) 



Indeed, the anomalous dimension is proportional to the Melhn transform of the gluon sphtting function 

CIsCf f 



7(cj) 



Pgg{z) z"^ dz . 



Energy conservation then imphes that 



7(1)=^^ rP,,iz)zdz = 0. 



(63) 



(64) 
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where 



lo = \l^-as (69) 



is the saddle point. Substituting (68) into te6p and integrating over the saddle point gives 



^ Q s 

The most important correction due to energy conservation requirement is steeper dependence of 
the scattering amplitude on r. 



2. Saturation momentum 

To determine the saturation momentum, we need to find a set of lines in the y, S, plane along 
which the amplitude is constant. In the DLA approximation this is equivalent to the requirement 



that the phase (68) be constant, i.e. 2y/a^y^ — ^(1 + as) = 0. Denoting solution to this equation 



as S,s{y) we obtain 

Ql = Qle^' = Ql eoTSTF . (71) 



iasy 



Energy dependence of the saturation momentum becomes more gradual compared to the LO. 

A more accurate evaluation of the saturation momentum requires solving the following two 
equations [59j: 

_ 2asx{l)y 



+ 7^ - ^ = (72a) 

2 +e = 0. (72b) 



1 + 2a,x(7) 

2asX'{l)y {2asfx{^)x'{l)y 



dj 1 + 2asX{l) (1 + 2asX(7)) 
The first one determines the line on y,^ plane where the amplitude is stationary, while the second 
one fixes the trajectory of the steepest descend [59] . Eliminating y and S, from these equations we 
end up with an equation for the saddle-point 7sp: 



Employing (rol) we write 



// N , 1 / ^ 2asX 7sp X (7sp) /»o^ 

X (7sp + z X 7sp) = -, , „- — J ^ • (73) 

1 - 7sp 1 + 2asX(7sp) 



x(7) = V'(l)-^V'(7)-^^(l-7), (74) 

x'(7) = -^V''(7) + ^^'(l-7). (75) 
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Saddle point in the LO is obtained as the solution to (73) in the a^ — )• limit. Hence, dropping 



the r.h.s. of (73) we obtain 7sp = 0.37. In the NLO approximation 7sp depends on a^ as shown in 
Fig. W[a). As ots increases 7sp decreases and becomes closer to the experimental data. For a given 



as ( 72 ) implies that 



Ql = Qlo exp 



1 2asx{lsv)y 
1 - 7sp 1 + 2qsx(7sp) 

x(7bp) 



q2 ^2asyh{as) 



(76) 



Particularly, at the LO h{as) = fi!Z' = 2-44 independently of as- In Fig. 2fb) we show the NLO 



behavior of h as given by (l76|) and its DLA given by ( 71 ). Again we observe that the NLO correction 



makes the energy dependence of the saturation scale more gradual. This is understandable because 
the energy conservation reduces the phase space available for gluon emission. 




(a) 



{b) 



FIG. 2: (a) Solution for the saddle point equation (73) 7sp(as): solid blue line is NLO (energy conservation), 
dashed red line is LO. (b) Function h{as) defined in ([76|: solid blue line is NLO, dotted (purple) is its DLA 



(71 1 and dashed (red) is LO. 



3. Saturation regime 



In the saturation region, ( 58 ) reads 

5iV(r,b,y) _ 

= a.o 1 



d 



dr 



12 



rl2 



{iV(r',b,y)-A^(r',b,y)A^(r,b,y)} 



(77) 



dy "° V dyj J2/QI r'- 

We expect that the scattering amplitude will approach its unitarity limit as y — )• 00. Therefore, 



we are looking for a solution to (77) in the form 

N = l-S 
where S* <C 1 is an element of the scattering-matrix of dipole r. Now 



(78) 



(79) 
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We are interested in the scaling solution viz. we are looking for a solution in the form S{r, y) 
S{T{r, y)) where 



2/-i2 



^O-sV 



T = Hr'Qt) = In(r^gfo) + n ^^y 



(80) 



and we used (71). Introducing a new parameter that determines rapidity dependence of the satu- 



ration scale (in the DLA) 



4a., 



:i+a.)2 



we write ( 79 ) as 



as; 

57 



(agAr — \) = as{T — X)S . 



(81) 



(82) 



It is easily integrated with the solution 

S{T) = Soei{l-a,T)^-\ (83) 

where Sq is an integration constant that is determined by matching with the solution in the dilute 



regime. This is similar to the solution derived in [64j- Note, that (83) is applicable only at 



1 < T < l/ttg. Solution (83) is exhibited in Fig. ^ 



1.0 
0.8 

^ 0.6 

I 

0.4 
0.2 
0.0 



FIG. 3: Solution to the LO (dashed red line) and the modified (solid blue) BK equations deeply in the 
saturation region 1 < r < l/oLg. The initial condition is S* = 0.9 at t = 1. 



B. Dipole density 



We proceed with the analysis of the NLO effects related to the energy conservation in the dipole 
density. Using the result of the Sec. VIA we obtain in place of ([9]): 



/(E,q,y) 



vrg 



du exp 



2asX{i')y \ rrqV^i^ V {\ - iv) 

1 + 2asx{y) J V 2 ; T{\+iv) {2iu + 1)2 



(84) 
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where v* satisfy 1 + 2asx(^*) = 0. Similarly to our discussion in Sec. Ill A, we would like to find 



asymptotic expressions for / in various kinematic regions. Since the integrand in (84) is a steeply 
falling function of u we can replace the limits of integration by v* = ztcc.^ 



1. asV ^ In ^. Expression in the exponent of (84) can be approximated as 



2asX{y)y ^ {af - l)y 14C(3)a,2/^2 



l + 2asx{^) 



a 



(0) 



a 



(0)12 



(85) 



We see that the pomeron intercept became ap = 2 — l/op , while the "diffusion constant" 
has increased by l/[ap ]^, i.e. growth of / with rapidity has slowed down, while diffusion has 
speeded up. The later observation has profound implications on diffractive gluon production 
(see [55H7T] for in-depth discussion). For Ug = 0.4 the intercept is Up =1.5 (compare with 



Op = 2.1), which is in better agreement with the data. Eq. (13) is modified as follows 



/(E,q,y) 



a 



(0) 



v(l) 



4°)]2ln2l| 



gVl4<(3)a,y 



("p -^)y Q 14C{3)as (y-H) 



(86) 



2. rq < 2 and In ^ >> a^y. Expanding x ^ 1/(1 — 2iz^) we find the saddle point at 



2iui = 1 + 2as 



l2asy 



In 



_2_ • 
rq 



Integration over the saddle-point and assuming In ^ ^ y/ds yields 



fii,<i,y) 



{rq/2) 



2as 



{2asyln^y' 



2. /2cJsj/ln 



1- ./2a ,Mn- 



'y rq 



(87) 



(88) 



3. rq > 2 and a^y <C In ^ <C y/os- Now, another pole in x dominates x ^ 1/(1 + 211^) with 
the result 

1 (2/rg)2"= 



f{i,<i,y) 



29'V^(2a,ylnf)'/" 



1- J2a4ln1 



^2^2as3/ln^ 



(89) 



Note, that in both cases (88) and (89) the momentum dependence of the leading twist is modified 
by an additional power 2as- This can have important consequences at high Q^ and/or kx- We are 



discussing this in more detail in Sec. VIII 



^ Note, that we keep u* finite for the purpose of the numerical integration in Sec. VIII 
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VII. NUCLEAR MODIFICATION FACTOR 



The nuclear modification factor is defined as 



Jd%^^^^ 



R^,^ = ^ ^^^g^ . (90) 



In the logarithmic approximation ( 56 ) implies that the cross section for inclusive gluon production 



in DIS on a heavy nucleus is simply proportional to the cross section for inclusive gluon production 



by dipole of size r = 2/Q. Consequently, the nuclear modification factor (90) can be approximated 
by 

In the same approximation, pA scattering can also be approximated as the qq + A one provided 
that we are interested in inclusive processes not too close in rapidity to the proton or nucleus 
fragmentation region [Til- Atomic weight A and rapidity y dependence of incluisve cross section 
in pA collisions at the leading logarithmic order was discussed in great detail in [34j and we refer 
the interested reader to that paper. Here we will focus on the logarithmic derivative of the nuclear 
modification factor defined as 

J =^2^. (92) 

R^*A dy 

Outside the saturation region this observable is proportional to the difference between the anoma- 
lous dimension of the gluon distribution in the nucleus 7 and the one in the proton 7^. If the 
coherence effects were negligible, the two anomalous dimensions would have been identical. This is 
not the case according to the theory of gluon saturation. As the result, the NMF is suppressed even 
at kx > Qs- Thus J is especially sensitive probe of the mechanism that leads to the suppression 
of the NMF for hadron production at small x. 



Let us relate J to the difference of anomalous dimensions 7 — 7^. It follows from (90) that 



/ - — 1 R - — 1 ^^^*^ _ A ^^^*^ 
dy ^ dy d'^kdy dy d'^kdy 



Using (91) and (30), (31) and assuming that the 6-dependence factors out we derive 



-lnR,,A^-lng |^^„ - - 



l"^7*A~7^1n5^|b=o-7^1n5P|^^p, (94) 



where g is the inclusive qq + A cross section modulo a constant factor, see Sec. IV We assigned 



superscripts A and p to g to indicate the two cases: A > 1 and A = 1 respectively. In the following 
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we will omit the specification that g is taken at zero impact parameter. Outside the saturation 



region we can employ the Mellin representation for Nq (23) and / i9v, substitute them into (31), 



take the LLA limit and obtain up to a pre-exponential factor 
g oc Qor exp 



2asx{i^o){Y -y) + 2iuo In ^ + 2a,x(/i^)y + 2i/i^ In ^*° 



rp 

y 



p 



(95) 



and analogously for g^. Here uq, /Iq are the saddle points in the Mellin transform of / and Nq 



respectively. The omitted pre- factor in (95) depends on momenta only logarithmically. Momentum 



p stands for either Q or k depending on the kinematic region of interest. It is straightforward to 
verify that g and g^ obey the equations 



dg^ 
dy 



2a.[x(M^)-x(i^o)]/, 



dgP^ 
dy 



2as[x{f^l) - xii^oW 



This is just the Mellin transform of the BFKL equation. Plugging (96) into (94) we derive 



J = 2a, [x{4) - xif^D] « 2a,x'(7o') U " To") • 



(96) 



(97) 



x'il) is given by (75) and the saddle point 7q satisfies (73). 

Consider a few examples. Denote p = uiax{kT,Q}- In the region In ^f- ^ a^y we have (see 



e.g. (15) and (24)) 



1 



1 



X 



1 - 2ifi 27 



with the saddle point 



2asy 



^^4(--^o)4J,,,,^ 



1 &A , 1^ 
2\ In^ r^ 21n£ 



(98) 



(99) 



A 



A 



7^ is obtained by setting Qsq = A. We see that in this kinematic region 7^ < 7 . By dint of (98) 
x'il) < implying that J < 0. More precisely, 



J 



-a.o 



m ^ 



2asy In f 



(100) 



In the saturation region In y^ ^ a^y, 7"^ effectively tends to zero as the dipole scattering amplitude 
saturates at unity. Therefore, in that region 7 < 7^, while x ^ 2d- 1 ■ Hence x'il) > implying 
that again J < 0. Finally, in the diffusion region x ~ 21n2 — 7C(3)z/^ and we similarly obtain 



J 



In £^ln 2^ 
m ^ m ^ 



(101) 



7C(3)a,y2 ' 

Negativity of J in all kinematic regions signifies the decrease of the inclusive cross section as a 
function of rapidity. The rate of the decrease depends on the absolute value of J. 
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FIG. 4: Comparison between the LO and NLO calculations of k^ — d'^kd '^ ^^ ^ function of kx at two values 
of coupling (a) cis — 0.3 and (b) as — 0.15. 



VIII. NUMERICAL ANALYSIS 



The numerical calculation of the inclusive hadron production is performed using 
Eqs. p^,(p|),([IT|),(10). We employed the bGCG model [67] for the forward dipole-nucleus scat- 
tering amplitude. The bCGC model is reviewed in Appendix. Function / is calculated using 



formula (84). The gluon spectrum is then convoluted with the LO pion fragmentation function J^q 



as follows 



day. 



d?kdy 



dz dac 



z^ d'^kdy 



ik/z)TGiz,k). 



(102) 



The fragmentation function is given in [80j. The total rapidity interval is taken to be y = 10, 



which is equivalent to a; = e' 



-Y 



4.5 -10 ^. The range of photon virtualities that we consider is 



Q^ = 2 — 37 GeV^. This kinematic region can be probed at the proposed Large Hadron electron 
Collider and its low Q^ part at the Electron Ion Collider [81]. The rapidity interval y from the 



nucleus to the produced gluon is related to xp, a variable used in differctive DIS, as xp = e~^. 
We consider y in a narrow interval 5 < y < 7 allowed by our formalism. At larger x and/or xp 
the validity of the leading logarithmic approximation that we employ becomes uncertain. 

The results of our calculations are shown in Figs. Q- ([8]). The NLO calculation shown in the 
figures refers to the part of the NLO terms that are responsible for energy conservation. In Fig. 4|5 
we plot the inclusive cross section normalized in the same way as the structure function 

dF2{x,Q^;y) ^ 1 Q^ da^*^ix,Q^;y) 
d^kdy aem 47r2 d'^kdy 

We observe that inclusive gluon production at NLO is suppressed compared with the LO case. 
This is because the anomalous dimension of dipole density at NLO is smaller compared with that 
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FIG. 5: Inclusive spectrum fc^ — d^kdy °^ (^) pio^is, (b) gluons as a function of /c^. 

of LO, as can be seen in Fig. [2] This is expected since energy conservation constrains the phase 
space available for hadron production. In Fig. [4] we demonstrate that the difference between the 
LO and NLO calculation is smaller at smaller values of coupling. 

We see in Fig. l5[b) that at small k^, the gluon production cross section follows l/k"^ behavior. 
Indeed, l/k"^ comes from the Lipatov vertex, whereas the gluon distribution in the nucleus is 



saturated and hence depends on momentum kj- only logarithmically. This is seen in (11) where at 
small kT the integral tends to a constant leaving the l/A:^ pre-factor in front. Modification of the 
gluon spectrum due to fragmentation can be inferred by comparing Fig. pla) and (b). 

The cross section grows with Q^ and xp logarithmically; both dependences are much steeper at 
the LO than in the NLO. We also note that energy conservation correction substantially reduces the 
cross section. However, the functional form of the A^T^-spectrum does not change in the kinematic 
region that we studied, as we checked explicitly. We attribute this to that fact that the dominant 
contribution to the Mellin transform stems from anomalous dimension 7 ~ 1/2 in both cases. We 
expect that at much larger Q and k^ the NLO A^j^-spectrum becomes steeper than those in LO due 
to additional factors l/Q^"" or l/kr^" . However, assumptions of our model restrict our calculation 
only to the semi-hard values of transverse momenta. 

The largest uncertainty in our numerical calculation of hadron spectrum comes from the over- 
simplified treatment of nuclei geometry. Instead of integrating with a realistic nuclear thickness 
T(6) we approximated the nuclear density by the step-function. Based on our previous experience 
with this type of numerical calculations we expect that a more accurate treatment of the nuclear 
density will only affect the overall normalization of the cross section. From this perspective the 
ratios of the inclusive spectra should not be much affected by this uncertainty. 

Our calculation of the Nuclear Modification Factor (NMF) as a function of kT for Au {A = 197) 
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FIG. 6: Inclusive hadron spectrum k^ — d^kdv ^® ^ function of (a) y, (b) Q'^. 

and Ca {A = 40) is displayed in Fig. u\ The general feature of NMF is suppression at low kj- and 
enhancement at larger Ut (the later is often referred to as the Cronin effect). This is in contrast 
with the hadron production in pA scattering where the Cronin effect gives way to the suppression 
at all /ct's provided that the hadron rapidity y is large enough. The reason for this difference is that 
whereas pA scattering can be approximated by dipole- nucleus scattering [7Tj, ^* A interaction is a 
superposition of many dipole- nucleus scatterings with different dipole sizes r, see ([3]). At small kx 
NMF for dipoles of all sizes is suppressed [34J and therefore we observe suppression of the resulting 
R for DIS. On the other hand, the fact that i? > 1 at large k^ implies that the inclusive cross 
section in that region is dominated by dipoles whose individual scattering on the nucleus exhibits 
Cronin enhancement, i.e. they are not much effected by the low-x evolution. Presence of such 
dipoles is ensured by evolution of the dipole density n, which happens \iY — y ^ 1. Comparing 
FigsJ7^(a)-(c) with (d) we note that due to fragmentation, NMF of hadrons is much slower function 
of Q^, y and kj- than NMF of gluons. Additonally, fragmentation shifts the value of the transverse 
momentum at which NMF crosses unity towards lower /ct- 

Another feature seen in Fig. u\ (especially (d)) is that suppression of NMF at low kx and its 
enhancement at high kx increases with the photon virtuality Q^ . To understand the Q^ dependence 
of the NMF we note that a typical term in its twist expansion looks like 

where n > 1 is an integer number. It implies that 

dR 



(104) 



-n{-iA - lp)R- 



(105) 



91nQ2 

At large kx 7"^ > 7^ thus q'(^q2 < 0, whereas at small kx J^ < 7^ thus qiJq2 > 0. This is indeed 
what we observe in Fig. [7} Dependence of NMF on y can be explained similarly. 
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FIG. 7: Nuclear Modification Factor as a function of fc^ for (a)-(c) hadrons at various A, y and Q'^; (d) 
gluons. All calculations include the NLO effects. 



Fig. ^ displays the logarithmic derivative of the NMF J defined in (92). As we argued in 
Sec. VII| this quantity is proportional to the difference between the anomalous dimensions of the 



gluon distribution function in nucleus and proton, see(97). Our analysis in (100), (101) indicates 



that J is negative and decreases as the hadron rapidity y increases, which is indeed seen in Fig. [8| 
Similar trend has been noticed m.pA collisions in [78j. We can also see the effect of fragmentation on 
J by comparing Fig.^a),(b) with (c),(d). It is interesting that fragmentation completely erases the 
kx dependence, while leaving the y dependence qualitatively similar. We think that experimental 
investigation of J is of great interest as it emphasizes the difference between the (linear) gluon 
evolution in a heavy nucleus and in proton. 

IX. SUMMARY 



In this paper we studied the inclusive hadron production in DIS scattering at small x using the 
dipole model [38j. We presented the analytical formulas for the cross section in various kinematic 
regions and discussed the role of the energy conservation, which is perhaps the most important NLO 
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FIG. 8: Logarithmic derivative of NMF for dipole-nuleus scattering as a function for kx for (a),(b) gluons, 
(c),(d) hadrons. dipole size r, total rapidity Y and nuclear wight A are indicated on each plot. All 
calculations include the NLO effects. 



correction. Employing the modified BK equation suggested in [64tl65|. we derived the corresponding 
correction to the pomeron intercept and found that it is numerically closer to the phenomenological 
value than the LO result. We also computed the high energy asymptotic of the forvi^ard dipole- 
nucleus scattering amplitude. 

Motivated by possible low x DIS experiments with heavy nuclei [8T] we performed numerical 
calculations of the DIS inclusive cross section using the bCGC model [67]. The results are shown 
in Figs. |4]-[8J We noticed that the NLO effects generally tend to reduce the cross section and make 
it weaker function of its arguments as compared to the LO result. The nuclear modification factor 
exhibits suppression at low kx and enhancement at higher kx even at the largest hadron rapidities 
that we can address in our approach. To understand dependence of the NMF on rapidity better we 
introduced the logarithmic derivative of NMF J and showed that it is proportional to the difference 
between the anomalous dimension of the gluon distribution function in nucleus and proton. Since 
this difference is non- vanishing only due to coherence effects, J provides a direct measure of the 
effect of coherence on inclusive cross section. Figs. 7|8 show dependence of NMF and J on the 
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photon virtuality Q^, x and hadron rapidity y. We believe that our results may be helpful for 
experimental investigation of the low-x regime of QCD in DIS. 

Acknowledgments 

This work was supported in part by the U.S. Department of Energy under Grant No. DE-FG02- 
87ER40371. 

Appendix A: bCGC model 

We performed the numerical calculations using the bCGC model of the forward dipole scattering 
amplitude [67J. We treat the nuclei and proton profiles as step-functions; the saturation scales are 
assumed to scale with A a.s Ql cc A^'^ . The advantage of this model - besides its compliance with 
the known analytical approximations to the BK equation |79j - is that its parameters are fitted to 
the low X DIS data. The explicit form of the scattering amplitude N is given by 

N{T,0,y)= { V 4 ; ' (Al) 

[ l-exp[-aln2(6rQ,)], rQ, >2, 

where Q^ is the the quark saturation scale related to the gluon saturation scale Q'^ - which we 
have called simply the 'saturation scale' throughout the paper - by Q^ = (4/9) Q^. Its functional 
form is 

Q2 = ^l/3^A gAj/ ^A/2 Qgy2 ^ (^2) 

where s is the square of the center-of-mass energy and y is rapidity with respect to the central 
rapidity. The anomalous dimension is 

^ = ^- + wi \ ^ N In (^) ■ (A3) 

cX{ln^/s + y) \rQsJ 



The gluon dipole scattering amplitude can be calculated using (18). Parameters 7^ = 0.628 and 
c = 9.9 follow from the BFKL dynamics [79j, while Mq = 0.7 and A = 0.28 are fitted to the DIS 
data. Constants a and b are uniquely fixed from by the requirement of continuity of the amplitude 
and its first derivative. 
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